Dynamical study on polaron formation in a metal/polymer/metal structure 
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By considering a metal/polymer/metal structure within a tight-binding one-dimensional model, 
we have investigated the polaron formation in the presence of an electric field. When a sufficient 
voltage bias is applied to one of the metal electrodes, an electron is injected into the polymer chain, 
then a self-trapped polaron is formed at a few hundreds of femtoseconds while it moves slowly under 
a weak electric field (not larger than 1.0 x 10 4 V/cm). At an electric field between 1.0 x 10 4 V/cm 
and 8.0 x 10 4 V/cm, the polaron is still formed, since the injected electron is bounded between the 
interface barriers for quite a long time. It is shown that the electric field applied at the polymer 
chain reduces effectively the potential barrier in the metal/polymer interface. 

PACS numbers: 71.38.-k; 72.80.Le; 73.40.Ns 
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I. INTRODUCTION 

Recent years, organic electronic devices (OEDs), e.g., 
light-emitting diodes (LEDs) and field-effect transistors 
(FETs), are attracting considerable interest because they 
have processing and performance advantages for low- 
cost and large-area applications^ In these devices, or- 
ganic polymers are used as the light-emitting and charge- 
transporting layers, in which the electron and/or hole are 
injected from the metal electrodes and transported un- 
der the influence of an external electric field. Due to 
the strong electron-lattice interactions, it is well known 
that additional electrons or holes in conjugated polymers 
will induce self-localized excitations, such as solitons 2 
(only in irans-polyacetylene) and polarons^* As a result, 
it has been generally accepted that the charge carriers 
in conjugated polymers are these excitations including 
both charge and lattice distortion^ The formation and 
transport of such charge carriers are believed to be of 
fundamental importance for these OEDs. 

There have been extensive studies on soliton and po- 
laron dynamics in conjugated polymers£i£iL& under the 
influence of external electric fields. It is shown that soli- 
tons as well as polarons keep their shape while moving 
along a chain. Solitons are shown to have a maximum 
velocity 2.7v 8 , where v s is the sound velocityA^ The 
situation will be different for polarons, which has been 
shown to be not created in electric fields over 6 x 10 4 
V/cm due to the charge moving faster and not allowing 
the distortion to occuri A recent study by Johansson 
and Stafstrom 8 deals with the polaron migration between 
neighboring polymer chains. The numerical results show 
that the polaron becomes totally delocalized, either be- 
fore or after the chain jump for the electric field over 
3 x 10 5 V/cm. There are also studies on the charge 
transport through metal/polymer interfaces, e.g., the re- 
sistance of organic molecular wires attached to metallic 
surface 10 , Schottky energy barriers in metal/organioii, 
and the dynamics of charge transport in a short oligomer 
sandwiched between two metal contacts^ 



In this paper, we present our results of a non-adiabatic 
dynamical study on the polaron formation starting from 
charge injection in a metal/polymer/metal structure in 
the presence of an external electric field. We would like 
to focus on charge and its induced lattice distortion in 
the polymer chain while the metal electrodes at the two 
ends serves as the source and drain for the extra charges 
within a tight-binding one-dimensional model. The elec- 
tron wavefunction is described by the time-dependent 
Schrodinger equation, in which the transition between 
instantaneous eigenstates are allowed, while the polymer 
lattice is treated classically by a Newtonian equation of 
motion. 6 

The paper is organized as follows. In the following sec- 
tion, we present a tight-binding one-dimensional model 
for the metal/polymer/metal structure and describe the 
dynamical evolution method. Our results are presented 
in Sec. Ill and the summary of this paper is given in Sec. 
IV. 



II. MODEL AND METHOD 

We consider a one-dimensional metal/polymer/metal 
structure that contains a polymer chain as well as two 
metal electrodes attached at its two ends. The Hamilto- 
nian consists of three parts, 



H = H e + H, 



latt 



The electronic part is 



H e = - ^t n (clc n+1 + h.c), 



(1) 



(2) 



t n being the hopping integral between site n and n + 1, 
that is, t n = tg in the two metal electrodes, t n — tg — 
a.{u n +\ — u n ), where u n is the monomer displacement 
of site n and a describes the electron-lattice coupling 
between neighboring sites in the polymer chain, as the 
Su-Schrieffer-Heeger (SSH) model^, and t n = t\ for the 



coupling between sites connecting the polymer chain and 
the metal electrodes. The polymer lattice is described by 



(a) 
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where K is the elastic constant and M the mass of a CH 
group^ The contribution from the external field is 



H ext =^V n (i)(ct Cn -l) 



(4) 



V„(t) being site-energies due to the applied voltage bias 
and electric field. At the left metal electrode a voltage 
bias is applied for the charge injection, V n (t) — V(t). 
At the polymer chain an uniform electric field E(t) is ap- 
plied, V n (t) — \e\E(t)[(n— n )a+u n ] with e being the elec- 
tron charge, no the first site of the polymer chain, and a 
the lattce constant. At the right metal electrode, a field- 
free area, the site-energies V n {t) are chosen as V ni (t), n\ 
is the last site of the polymer chain. The spin index in 
the electron operators is omitted since the electron in- 
teraction will not be considered here. A similar static 
model is used very recently to investigate the ground- 
state properties of ferromagnetic metal/conjugated poly- 
mer interfaces^ 

We consider a finite system containing a 200-monomer 
polymer chain and two 100-site metal electrodes at- 
tached at the two ends. The model parameters are those 
generally chosen for irans-polyacetylenei^ to = 2.5eV, 
a = 4.1eV/A, K = 2IeV/A 2 , M = 1349.I4eVfs 2 /A 2 , 
and a — 1.22 A. The coupling between the polymer and 
the metal is set to t\ — l.OeV. Before we go further for the 
dynamical evolution, we determine the static structure of 
energy levels in a constant external field (V(t) = Vo and 
E(t) = Eq) in this section. 

The total energy is obtained by the expectation value 
of the Hamiltonian (1) at the half-filled ground state \g), 

K 

E t = (g\H e + H ext \g) + — 22(u n+1 - u n ) 2 . (5) 

n 

The electronic states are determined by the electronic 
part of the Hamiltonian (1) and the lattice configuration 
of the polymer {u n } is determined by the minimization 
of the total energy in the above expression, 



-'n+l 



-^KPn,n+l + Pn+X,n) H f^Pn.n + A, (6) 

K K 



where A is a Lagrangian multiplier to guarantee the poly- 
mer chain length unchanged, i.e., X)n( u «+i ~ M «) = ®- 
p n ,n' is the element of density matrix, which will be given 
below. 

Due to the coupling between the polymer and the 
metal electrodes, the electronic states in the polymer 
and the metal electrodes are mixed. By calculating 

the wave function possibilities Tjj^\= Xme« I0p( n )| 2 ) f° r 
each state /x, where ^(n) is the wave function at site n 
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FIG. 1: Schematic energy level diagrams of the 
metal/polymer/metal structure in the absence (a) and 
presence (b) of external fields. E g is the band gap. The two 
dashed lines indicate the localized levels if a polaron exists 
in the polymer chain. Vo is the voltage bias applied at the 
left metal electrode and Eo is the electric field applied at the 
polymer chain. 



and k is a set of sites (the left metal electrode, the poly- 
mer chain, or the right metal electrode) in the system, we 
can find that each electronic state can be clearly distin- 
guished to be in the polymer or the metal electrodes, ft 
indicates that the coupling (ti/to) in the metal/polymer 
interfaces is small. 

In Fig. f , we show the schematic energy level diagrams. 
In the absence of external fields, the polymer chain is 
dimerized for a half-filled system, the energy band, as in 
a single polymer chain, has a gap of about E g = 1.32eV, 
and the metal electrodes have quasi-continuum bands 
with the Fermi level almost at the middle of the gap. The 
dashed lines indicates the localized levels if a polaron ex- 
ists in the polymer chain. E^ « ±0.45eV, indicates the 
binding energy of a polaron to be about 0.2IeV. It has 
been known that the energy gap of MEH-PPV is 2.4eVi^, 
a much larger value than the energy gap E g we obtained 
above, but the polaron binding energy is similar to that 
of MEH-PPV. 

In the presence of a constant external field (the voltage 
bias Vq and the electric field Eq), the energy levels will be 
changed as shown in Fig. 1(b). The voltage bias Vo ap- 
plied at the left metal electrode raises its energy levels by 
a value of about Vb while the energy levels in other parts 
are almost unchanged at Eq = 0. It has been showni^ 
that when the highest occupied electronic state of the left 
metal electrode is close to the bottom of the conduction 
band due to the applied voltage bias (Vb > 0.56eV), the 
system will be unstable and charges start to be injected 
into the polymer through the metal/polymer interface. 
It will be seen that the electric field E applied at the 
polymer chain favors the charge injection as a result of 
the incline of the site-energy at polymer chain. 

The voltage bias Vo could also serve as the metal work 
function. At Vo = 0, the energy difference between the 
electron polaron level in the polymer chain and the high- 
est occupied level of the metal electrode is about 0.45eV. 
For the contact with MEH-PPV and aluminum, the en- 
ergy difference is about I.4eV, and the value is about 
for that with MEH-PPV and calciumii So the value 
Vb = 0.57eV, we take in the following calculation, could 
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be thought to correspond to that for the Ca/MEH-PPV 
contact with the voltage bias about 0.2eV applied at the 
metal calcium. 

Now, we describe the non-adiabatic dynamical method 
that has been used for the dynamics of soliton-i and 
polaroni^ in an electron-lattice interacting system. The 
evolution of the electron wavefunctions depends on the 
time-dependent Schrodinger equation 

where the site index n runs over the whole chain. The 
lattice displacements are determined classically by the 
following Newtonian equations of motion 

Mu n {t) = K[u n+l {t)+u n _ 1 {t)-2u n {t)] 
+2a [p n) „ + i (t) - p n -x t „(t)] 
+\e\E(t)[p n , n (t) - 1], (8) 

where n runs only in the polymer sites. p nn ' is the ele- 
ment of the density matrix defined as 

Pn,n'(t) = ^4>*n^{t)fli<t>n', ll {t), (9) 

where /„ is the time-independent distribution function 
determined by initial occupation (being 0, 1, or 2). The 
coupled differential equations (7) and (8) can be solved 
numerically by use of the same technique in Refi£*&. The 
time step is chosen to be as small as 0.1 fs to avoid numer- 
ical errors. A small damping is also introduced to smear 
out the lattice vibrations, and the effect of the vibrations 
will be elaborated elsewhere. 

In the real calculation, we choose the external field 
to be turned on smoothly, that is, we let V(t) = 
V exp[-(t - t c ) 2 /tl] for < t < t c and V(t) = V for 
t > t c with t c being a smooth turn-on period and t w 
the width. As the voltage bias V(t), the applied elec- 
tric field E{t) is also turned on smoothly, i.e., E(t) — 
E exp[-(t - t c ) 2 /t 2 J for < t < t c and E(t) = E for 
t > t c with the same t c and t w . In calculations, we take 
t c = 30/s, t w = 25/s, and various values of voltage bias 
Vo and electric field Eq. 

III. RESULTS 

In this section, we present our results on the polaron 
formation starting from the charge injection from the 
metal electrode in the presence of an external electric 
field. For that, we put one extra electron on the low- 
est unoccupied state of the left metal electrode while the 
voltage bias and electric field are turned on. In the fol- 
lowing, we will focus on the evolution of the charge dis- 
tribution p n [= p n>n — 1] and the staggered lattice con- 
figuration y n [= (— l) n {u n — u n ~i)]- At the initial state 
(t = 0), the polymer chain has a dimerized lattice of 
y n ~ 0.08A but a few bonds at its two ends and p n = 
for all sites in the polymer chain. 
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FIG. 2: The charge distribution p n (solid circles) and the 
staggered lattice configuration y n (open circles) of the poly- 
mer chain at different times under a weak electric field Eq = 
5 x 10 3 V/cm. Vo = 0.57eV. 



First of all, we shows the charge distribution p n and 
the staggered lattice configuration y n (in unit of A) of the 
polymer chain at a few typical times under a weak electric 
field E — 5 x 10 3 V/cm in Fig. 2. The voltage bias ap- 
plied at the left metal electrode is taken as Vo = 0.57eV, 
which is just above the minimum value of the bias for 
the charge injection in the absence of the electric fields 
For a smaller value of Vo, the charge will not be injected 
into the polymer since the energy of the highest occu- 
pied electronic level at the left metal electrode (Eho) 
is much lower than the bottom of the conduction band 
of the polymer (E g /2). For a larger value of Vo, more 
electrons will be transferred into the polymer since more 
occupied levels of the left metal electrode will be higher 
than E g /2. 

As the voltage bias applied at the left metal electrode 
increases to Vo = 0.57eV smoothly, Eho goes up to- 
ward the bottom of the conduction band of polymer. At 
t = 30fs, Eho reaches its maximum value 0.65eV, which 
is very close to the bottom of the conduction band. As a 
consequence, the coupling between energy levels mainly 
at the left metal electrode and at the polymer becomes 
stronger and stronger, the extra electron initially on the 
left metal electrode begins to be injected into the polymer 
gradually. From the Fig. 3, which shows the evolution of 
the total charges in the polymer chain and the metal elec- 
trodes, we can see that the charge injection takes a few 
hundreds of femtoseconds, which seems to be insensitive 
to the strength of electric fields. 

Furthermore, from Fig. 3(a), one can find that the 
charge does not enter the right metal electrode, which 
shows the electric field is too weak to overcome the bar- 
rier at the right polymer/metal interface, though we have 
made the calculation up to a few ms. The situation will 
be different for the electric field Eq > 1.0 x 10 4 V/cm, 
see Fig. 3(b), which is strong enough for the charge to 
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FIG. 3: Evolution of the total charges in the polymer chain 
(line 1) and the metal segments (line 2 and 3) at the electric 
field E = 5 X 10 3 V/cm (a) and E = 3.0 x 10 4 V/cm (b). 
V = 0.57eV. 



overcome the interface barrier and finally enter the right 
metal electrode in a few thousands fs. 

In contrast with the case of Eo = 0^ the polaron is 
formed not at the left side of the polymer but at the 
right side, as shown in Fig. 2, due to the action of the 
electric field. It will be seen more clearly from Fig. 4, 
which shows the evolution of the charge center x c in the 
polymer, which is defined as 



200- 



<D 
SI 

o 



150 



100 



1 . 


11,1111, 
E = 




- 


2. 


E = 5x10 3 V/cm 






13. 


E = 3x10 4 V/cm Ju, 






" 4. 


E = 7x1 4 V/ctjK^ 










2 




'\J vir 1 ^ 


1 












. i i 







250 



500 
time (fs) 



750 



1000 



FIG. 4: Evolution of the charge center x c in the polymer chain 
for various electric fields. Vo = 0.57eV. 
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FIG. 5: Same as in Fig. 2 but with E = 1.0 x 10 5 V/cm and 
Vo = 0.57eV. The inset shows the total charges in the polymer 
chain (solid line) and the metal electrodes (dash and dot-dash 
lines). 



Xc = ^2 n Pn/^2Pn, (10) 

n n 

where the summation is over the polymer sites. It can 
be seen that there are two regions, one is the charge in- 
jection at t = — 200fs, and the other is the polaron 
formation at t = 200 — 500fs. In the first region, the 
charge is injected into the extended states of the poly- 
mer and moves under the electric field. Since the local- 
ized lattice distortion has not yet formed in this region, 
the injected electron behaves as a free conduction elec- 
tron. As a comparison, we also show the curve for Eo = 
in Fig. 4, where the charge moves away from the inter- 
face just due to the energetic favorite. At a weak electric 
field, such as Eo — 5.0 x 10 3 V/cm, the charge center 
moves gradually to the middle of the polymer chain at 
t = 200fs, which corresponds with the charge being al- 
most uniform in the polymer. However, at an electric 
field with the strength Eo between 1.0 x 10 4 V/cm and 
8.0 x 10 4 V/cm, the charge will be spread quickly and 
then bounced back at the right polymer/metal interface, 



which is also the reason that the polaron will be formed 
at a position closer to the left interface under a stronger 
electric field. In the second region t = 200 — 500fs, the 
lattice distortion makes the polaron form with the charge 
being trapped inside while it moves toward to the right 
interface under the electric field. 

What we have found also implies that the polaron will 
not form in electric fields > 10 4 V/cm if the charge does 
not be bounded by the metal/polymer interfaces. This 
critical strength of applied electric fields is smaller than 
that obtained by Rakhmanova and Conwelfi who gave 
the value of 6 x 10 4 V/cm for the polaron formation but 
with the initial state of the electron being injected into 
the polymer chain with the lowest possible energy. 

In the case of the electric fields Eo > 8.0 x 10 4 V/cm, 
there will be more electrons being injected into the poly- 
mer as a result of the incline of the site-energy at polymer 
chain. In Fig. 5, we show the evolution of the polaron for- 
mation for Eo = 1.0 x 10 5 V/cm. It can be seen from the 
inset of Fig. 5 that the total charge at the polymer has 
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FIG. 6: Same as in Fig. 2 but with E = 2.0 x 10 5 V/cm and 
Vo = 0.3eV. The inset shows the total charges in the polymer 
chain (solid line) and the metal electrodes (dash and dot-dash 
lines). 



reached 2e at t — 250fs, when the lattice displays a few 
distorted wells at random locations. When the charge 
continues to be injected into the polymer at the period 
t = 250 — lOOOfs, the lattice forms a single deep well close 
to the right interface and the charge also starts to enter 
the right metal electrode from the polymer under the 
strong electric field. At around t — 1200fs, the charge in 
the polymer decreases to about one electron and a well- 
shaped polaron locates at the position close to the right 
interface. It lasts for a few thousands femtoseconds, de- 
pending on the strength of the applied field, before the 
charge enters completely the right metal electrode and 
the polymer gets back to the dimcrized state. 

The result for the electric fields E > 8.0 x 10 4 V/cm 
indicates that the electric field applied at the polymer 
chain reduces effectively the potential barrier at the 
metal/polymer interface. To show that, we consider 
the case of a smaller voltage bias applied at the left 
metal electrode. For example, we take Vq = 0.3eV, it 



is found that the charge can be injected into the poly- 
mer for the electric fields > 1.7 x 10 5 V/cm. In Fig. 6, 
we show the result of the calculation at the electric field 
E = 2.0 x 10 5 V/cm and the voltage bias V = 0.3eV. 
The physical process is similar but the polaron lasts only 
a short period since the charge enters the right metal 
electrode sooner under a stronger electric field, which also 
indicates that the potential barrier at the polymer/metal 
interface is effectively reduced. 



IV. SUMMARY 

In summary, we have investigated the polaron forma- 
tion starting from the charge injection by using a non- 
adiabatic dynamic method based on the time-dependent 
Schrodinger equation for the electronic wavefunctions 
combining the Newtonian equation of motion for the 
polymer monomer displacements. It is shown that the 
polaron is formed while it moves slowly only under a 
weak electric field Eq < 1.0 x 10 4 V/cm. But for a 
stronger electric field, the polaron can be still formed 
because the charge is not easy to be jumped away from 
the chain. This result implies that the low mobility of 
the charge transport in polymers is mainly due to the 
obstacle for carriers to travel through the polymer/metal 
interface and/or jump over a neighboring chain. Our re- 
sults also show that the electric field applied at the poly- 
mer chain reduces effectively the potential barrier in the 
metal/polymer interface. 
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